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1.0 INTRODUCTION 

TRW has conducted design and development of Precision Attitude Deter- 
mination System (PADS) technology for NASA Goddard Space Flight Center 
(GSFC) under contract NAS5-21111. This report documents the feasibility 
evaluation of an evolutionary development for use of a single-axis gimbal 
star tracker from prior two-axis gimbal star tracker based system applica- 
tions. Detailed evaluation of the star tracker gimbal encoder is also 
addressed. 

The report provides a brief system description, including the aspects 
of tracker evolution and encoder evaluation. System analysis includes eval- 
uation of star availability and mounting constraints for the geosynchronous 
orbit application, and a covariance simulation analysis to evaluate perfor- 
mance potential. Star availability and covariance analysis digital computer 
programs are included as Appendices. 


2.0 SYSTEM DESIGN 

The Precision Attitude Determination System (PADS) employs star 
trackers and gyros to determine inertial -referenced attitude to arc-second 
level accuracy. PADS system design employing a two-axis gimbal star tracker- 
has been previously documented [1]. Key system design differences in em- 
ploying a single-gimbal star tracker involve both hardware (e.g., tracker 
and associated electronics) and software. These differences and potential 
system design and performance tradeoffs are discussed in this section. 

2.1 System Configurati on/Def ini ti on 

The PADS star tracker, electronics, gyros, and computer/software shown 
schematically in Figure 2-1, are functionally and operationally similar to 
previous design [i]. The single-gimballed Star Tracker Assembly (STA) is 
used to track stars and provide a measure of the star line-of-sight (LOS) 
angle. The Sensor Electronics Assembly (SEA) provides the gimbal servo/ 
drive electronics and gimbal readout processing. The Gyro Reference Assem- 
bly (GRA) incorporates a configuration of strapdown pulse rebalanced rate 
integrating gyros plus electronics. A Reference Block Assembly (RBA) pro- 
vides a stable thermo-mechanical interface with the spacecraft to which the 
sensors are mounted and aligned relative to each other. The Digital Pro- 
cessor Assembly (DPA) comprises the computer - I/O to implement system soft- 
ware. 

The gyros provide high-bandwidth data which is compensated to account 
for systematic gyro errors and used to derive attitude. Attitude and gyro 
drift (bias) corrections are periodically estimated via star tracker mea- 
surements processed in a Kalman filter. 

The gimbal led tracker is operated to intercept defined stars of oppor- 
tunity. The gimbal freedom is such as to provide a single axis of mechani- 
cal scan. The tracker is mounted in the spacecraft such that the FOV moves 
in a plane generally defined as normal to the orbit plane, with gimbal free- 
dom chosen to assure adequate star sightings. Based upon a series of such 
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PADS Functional Block Diagram 
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sightings, the attitude estimate will converge to an accuracy consistent 
with system alignment and bias uncertainties; typically, steady-state is 
reached within 10 sightings or less. During this period, the system will 
reach thermal equilibrium. In-flight calibration may be utilized to deter- 
mine compensation for systematic errors. Observables include gyro input 
axis misalignments, scale factor uncertainties, and tracker misalignments., 

The configurated PADS system physical characteristics are summarized 
in Table 2-1, The single-axis gimballed tracker utilizes the identical Star 
Sensor Unit (SSU) and one-axis of the gimbal suspension/drive assembly of 
the two-axis tracker design [1], Likewise, the Sensor Electronics Assembly 
is of identical design, although it again provides drive and readout elec- 
tronics for only a single channel. Revised size, weight, and power charac- 
teristics have been derived accordingly, 

2.2 Design Evaluation 

Prior PADS star tracker development, although oriented toward a two- 
axis gimbal design, was limited in hardware development to a single-axis 
engineering model gimbal. This single-axis design is assumed used directly 
for the current PADS application. During the PADS star tracker development 
and test [23» however, problems with the encoding electronics were noted but 
not fully defined. In the following discussion, these problems are identi- 
fied and design modifications recommended. 

In testing the PADS system, certain anomalous behavior was noted in the 
breadboard Inductosyn angle encoding. Investigation was confined to the 
single speed channel, which was modified to read out all 14 available bits. 
There were three basic forms of anomaly as shown by the error curve in 
Figure 2-2. 

(1) A wide-range cyclic error which appeared to be in excess of 2° 
peak-to-peak over the available 90° of gimbal travel. 
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Table 2-1. Summary of PADS Physical Characteristics 



Wt (lb) 

Power (w.) 

Envelope (in.) 

Star Tracker Assembly (STA)* 

26 

7 

9.6 X 20.0 X 

Sensor Electronics Assembly (SEA) 

4 

8 

6x8x4 

Gyro Reference Assembly (GRA)^^^ 

9 

10 

6x6x5 

Reference Block Assembly (RBA) 

10 

— 

— 

Digital Processor Assembly (DPA)^^^ 

17 

27 

11.5 X 10.5 X 6.5 

Integration Hardware 

4 

— 

— 

Total System 

70 lb 

52 w 

21" X 26" X 20"^^^ 

* Single-axis gimbal design 

(1) Triad, Unheated 
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(2) Typical 

(3) Approximate total envelope occupied 
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(2) A fine-scale cyclic spatial modulation of the error function, 
with amplitude approximately 0.1® and periodicity about 1.4®. 

(3) Noise in the output data of about 4 bits peak-to-peak. 

The investigation of these problems was initiated, with emphasis on 
diagnosis. The system was set up, including the SEA breadboard, the gimbal, 
and the test set. An optical cube mounted on the Star Tracker case permit- 
ted the use of a T3 autocollimator. A rotary table enabled rotation of the 
gimbal base with a resolution in the region of a fraction of a minute, the 
gimbal being manually counter-rotated to the Autocollimator null for each 
table setting. Note that because the original null is relatively arbitrary, 
accuracy is here to be understood in the sense of linearity. Absolute ac- 
curacy (bias) depends on the initial Inductosyn mechanical alignment in ro- 
tation, and on total system stability thereafter. As a first step, the test 
findings were essentially verified in all three respects: a gross non- 

linearity, a fine-scale linearity ripple, and excessive output data jitter 
of about 16 counts. 

The spatial ripple error was suspected to be due to stray coupling from 
the 256-speed channel because of its periodicity. The amplitude-to-phase 
breadboard was a poor layout, with a considerable amount of physical inter- 
leaving of the single-speed and multi -speed elements. To check this, all 
the integrated circuits of the multi -speed channel were removed, although 
the preamps were left operating. The ripple error was reduced from about 
0.1® peak-to-peak to a value indiscernable under existing conditions. 

The wide-range error was tracked to a miswired capacitor in the 
amplitude-to-phase converter. Correction of this yielded a linearity error 
of approximately 0.5 degree peak-to-peak (Figure 2-3). This is similar to 
the Farrand data for the single-speed Inductosyn, although the shape of the 
developed error curve cannot be matched with that for the Inductosyn because 
of unavailability of Farrand 's definitions of electrical zero, positive 
angular rotation, and positive error. 



























































































In an attempt to characterize the noise behavior, random data samples 
at a given gimbal angle were recorded and plotted. A typical histogram of 
100 samples is shown in Figure 2-4. The horizontal quantization is one 
least-significant count, of an electrical revolution. In the 256-speed 
channel, a quanta would be 0.309 arc-seconds. 

The analog processing chain, consisting of the Preamplifiers, the 
Buffers, and the Ampl itude-to-Phase Converters appears to be generating AM 
noise. This chain has a relatively high gain and uses yA715's, which are 
not notable as low-noise devices. This AM noise in the sinusoidal phase 
signals would be mapped into phase noise by the Zero-Crossing Detectors. A 
scope examination of these signals did liideed show some 200 to 300 ns of 
phase jitter of the logic transitions of the lag signal with respect to the 
lead. 


To check this, a passive RC phase-shifting network was constructed and 
excited directly with the excitation voltage, so as to give two sinusoidal 
outputs, one slightly lagging the other. These were fed directly to the 
Zero-Crossing Detectors, bypassing the analog front-end elements. The phases 
were now extremely stable, but the encoded output data was still quite noisy. 

To obtain a noise comparison, the output data was D-A converted and the 
result displayed on a scope and photographed with a slow sweep. Figure 2-5 
illustrates the resultant noise indication. Each dot represents a successive 
data sample. Each row of dots is associated with a distinct data quantiza- 
tion level of one least significant (2'“^'^) unit. 

The data for the "noiseless" phase-shifter appears to have a width of 
16 quanta and a relatively sharp-peaked distribution. That for the normal 
system has width about 20 quanta wide and a somewhat broader distribution 
peak. 
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The conclusion must be that there is a second noise source associated 

with the encoding process itself. That was tracked to the VCO within the 

14 

Excitation Section. In order to obtain 2 resolution with a moderately 
high kHz) excitation frequency and a manageable clock frequency (2.8 
mHz), phase determination is done over 15 excitation cycles. The excitation 
is generated as 15/2048 of the clock frequency. This has the effect of lin- 
ing up the clock pulses in successive excitation cycles, as illustrated in 
Figure 2-6. Each excitation cycle contains 273 1/15 clock pulses, leading 
to the 1/15 precession shown. 

A phaselock loop is used to obtain the integer 15 portion of the 
excitation-to-clock frequency relationship. Considerable short-term phase 
jitter was observed in the VCO output. The phenomenon of output data noise 
due to this can be related as follows. Consider the vertical lines in Fig- 
ure 2-6 to represent a time gate which admits clock pulses to a counter. 

As shown, the gate opens and closes at precisely the same phase each excita- 
tion cycle. A count of one or zero might be registered depending on the 
arbitrary choice of mean phase. 

Widening the gate with the same mean phase would represent the result 
of Inductosyn rotation. A random phase jitter of both gate edges, with re- 
spect to one another from cycle to cycle would represent the result of the 
phase noise previously described. There would be a corresponding variation 
in the count from one block of 15 cycles to another. 

Now assume that the gate width remains precisely constant, as would be 
the case for the passive network, but that its mean position jittered ran- 
domly from cycle-to-cycle. It is likely that at times it could sweep along 
at just the proper rate to gather up 15 clock pulses. At other times, it 
might intersect none. At still other times it may sweep up scime inter- 
mediate number. 
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Clock Pulse 



Figure 2-6. Clock Phases in Successive Excitation Cycles 


Hence, with a completely "noiseless" analog processing front end, but 
a phase-jittering excitation, noise is still possible. This phenomenon is 
a hazard peculiar only to multi-cycle encoding, since a single-cycle encoder 
would not use a phaselock loop. 

Basically, the PADS encoding technique appears to have been rather suc- 
cessful. During design, it was necessary to make certain practical compro- 
mises in the interests of working out and demonstrating the basic concepts. 
Primary stress was placed on accuracy. This is still not assessed because 
of our present inability to subtract the Inductosyn errors. But even if 
one assigns the whole ±0.2 degree linearity error of Figure 2-3 to the 
electronics alone, this is equivalent to ±3 arc-seconds at the multi -speed 
level. This could probably be improved to the ±1 arc-second vicinity. 

Little attention was paid to circuit shielding in PADS, partly because 
of its expense, but mostly because of the desirability of exposed circuitry 
in a developmental effort. It seems that a packaging style using a number 
of enclosed shielding compartments will be necessary for the analog front- 
end elements, at least for the multi-speed channels. 

Noise performance tended to be ignored in the PADS design. The choice 
of the jjA 715 for the front end elements was dictated by a desired for a wide 
bandwidth. Gain precision, and hence large excess gain at 10 kHz, is essen- 
tial to accuracy. This device is not characterized for noise. Newer ampli- 
fiers are available with reasonable bandwidths and noise characterizations. 
Although they do not have the amount of open-loop 10 kHz gain, it is possi- 
ble that a more judicious gain distribution can be used. 

The requirement for a very quite phaselock loop was unsuspected, partly 
because of the subtlety of the noise phenomenon previously described. Future 
design will probably require a crystal VCO. 

The present noise performance of the PADS gimbal readout appears to be 
of the order of ±8 to 10 quanta, corresponding to ±2.5 to 3 arc-seconds, 
multi -speed. Assuming that this is composed of approximately equal parts of 
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analog front-end noise and VCO noise, that the front-end noise could be re- 
duced by 5 using quieter amplifiers, and that the VCO noise could be elimi- 
nated by using a crystal VCO, it is likely that the noise performance could 
be reduced to the order of ±1.4 quanta, a fraction of an arc-second. Note, 
however, that this estimate is made without consideration for those sorts 
of non-linear sampling effects which may occur when noise and quantization 
approach each other. 

Evolution of the PADS software for the single-gimbal tracker system 
is limited to those elements which relate to models, measurement, and con- 
trol of the star tracker. All other software remains unchanged. Thus the 
only elements effected are the catalog processing (star sort, selection), 
computation of the measurement residual, and derivation of the measurement 
matrix. The catalog processing reflects logic to acquire and track stars 
of opportunity (rather than at periodic intervals), and the measurement re- 
lated algorithms previously derived are, in this case, constrained by the 
assumption of zero outer gimbal angle. 
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3.0 SYSTEM ANALYSIS 


3.1 Star Availability Analysis 

Assessment was made of the tradeoffs in geometry and mounting constraints 
for use of the existing single-gimbal tracker in geosynchronous orbit appli- 
cation. The key tracker constraints that must be satisfied include: 

• Star magnitudes brighter than m^ = +3.5 

• Tracker gimbal freedom ^ ±60 degrees 

• Sun avoidance angle of 45 degrees 

From a spacecraft systems application standpoint, certain other constraints 
and/or design guidelines were established, viz. 

• Only one single-gimbal tracker need be active 

• Minimize FOV geometry interface requirements 

• Configuration well suited to operation during all- times of year 

A star availability digital computer program was developed which orders 
stars (by orbit angle) in sequence of observation as constrained by the 
gimbal FOV, The orbit angle between successive star sightings becomes a 
measure of star availability. The digital program (Appendix A) includes 
specification of orbit parameters, tracker performance and mounting para- 
meters, and input from a separate master star catalog file. Table 3-1 sum- 
marizes the results of a series of computer runs. 

From the tabulated results, it appears that the best mounting geometry 
is for the nominal boresight (zero-gimbal angle) to be directed at one of 
the two poles. This minimizes the effects of sun interference (greatest for 
boresight N in summer, for boresight S in winter). Some flexibility in star 
mounting and selection appears achievable. Although choice of ±45° gimbal 
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Table 3-1. Single-Axis Giinbal Tracker Star Availability 
(Geosynchronous Orbit) 


Boresight 

Mounting 

Star 

Magnitude 

(minimum) 

Gimbal 

Freedom 

(deg) 

Time of 
Year 

Evaluated 

Results/Comments 

Zenith 

3.5 



Significant sun interference problem, particularly near equinoxes 

45®N 

3.5 

30 

Spring/ 

Sumer 

Spring - 69 stars, gap of 40° (about sun) 
Summer - 53 stars, gap of 119° (about sun) 

45°N 

3.5 

45 

Spri ng/ 
Summer 

Spring - 86 stars, gap of 40° 
Summer - 67 stars, gap of 106° 

N 

3.5 

45 

Spri ng/ 
Summer 

Spring - 58 stars, maximum gap of 27° 
Summer - 48 stars, maximum gap of 27° 

H 

3.5 

60 

Summer 

74 stars, few scattered gaps of 12-15° 

U 

4.0 

45 

Summer 

78 stars, scattered gaps of 10-12° 

45^5 

3.5 

30 

Winter 

80 stars, maximum gap of 75° (about sun), other scattered gaps of 10-22° 

45*S 

3.5 

1 

.. .. .j 

45 

I V/i nter 

' 

96 stars, maximum gap of 71° (about sun), other scattered gaps of 10-22° 

S 

3.5 

30 

All 

40 stars, scattered gaps of 10-23° 

S 

3.5 

45 

Spring 

104 stars, scattered gaps of 10-17° 




freedom appears reasonable, opening the freedom to ±60° improves performance 
potential whereas reduction to ±30°, while potentially degrading performance 
somewhat, minimizes spacecraft FOV constraints. It was also noted that se- 
lection of slightly dimmer stairs (e.g., m^ = +4.0) did, in fact, somewhat 
increase performance potential. 

Figure 3-1 is a histogram which summarizes tradeoffs for cases of bore- 
sight South. A baseline system configuration is selected with the boresight 
South and gimbal freedom of ±45 degrees. In geosynchronous orbit, 15 degrees 
of orbit motion corresponds to an hour. Thus, the worst case time between 
stars for the given orbit will be somewhat greater than one hour. 

3.2 Error Analysis 

PADS using a single-axis gimbal tracker will have performance character- 
istics which differ somewhat from a system using a two-axis gimbal tracker. 
This comes about primarily from the imposed limitation to utilize stars of 
opportunity rather than being able to seek stars at periodic (arbitrary) 
intervals. As a result, performance becomes more dependent upon gyro be- 
havior, since updates will not be available as often (as noted through the 
star availability analysis of the preceding section). TRW, in recognizing 
this, has done considerable work in the modelling and analysis of gyro ref- 
erence systems [ 3 ], results of which may be applied directly. This section 
summarizes the system and component errors - the system error budget is 
outlined in Table 3-2. The star tracker errors and gyro data for drift and 
bias uncertainty are based upon data and models developed in the sections 
which follow. Gyro alignment and scale factor errors assume the full un- 
certainty is realized with a uniform rate of change over 20 minutes. The 
average update period of 20 minutes implicit in the results for gyro error 
is consistent with the star availability in the preceding section. The RSS 
of these error sources betters the goal of 0.001 degree (= 3.6 arc-second), 
la, per axis. 


ZENITH 



SUMMARY OF REPRESENTATIVE CASES^'gEOSYNCHRONOUS ORBIT, ZERO 
BORESIGHT SOUTH (ALONG ORBIT NORMAL) 



^ data available for many cases. INDIVIDUAL stars and PERIOD 
BETWEEN SIGHTING IDENTI FI 


Figure 3-1 



Table 3-2. System Error Budget 
(Single-Gimbal Tracker) 


ERROR SOURCE/BUDGET 


CONTRIBUTION (la, per axis) 


STAR TRACKER 


BIAS 


0.44 sec (1 a) 


0. 44 sec 


RANDOM , 1.31 sec (I 0 ) 


1.31 sec 


NOISE 


1.2 sec 


1.2 sec 


SUBTOTAL, RSS 


1.83 sec 


GYRO REFERENCE 


RANDOM DR-IFT 


BIAS UNCERTAINTY 


<2 sec over 20 Minutes 


0. 001°/Hr 


1.2 sec 


ALIGNMENT STABILITY 


1 0 sec 


0.45 sec 


SCALE FACTOR STABILITY , 100 ppm 


0.9 sec 


PULSE QUANTIZATION 


0. 01 sec 


SUBTOTAL, RSS 


2. 54 sec 


COMPUTATION 


1.0 sec 


TOTAL, RSS 


3.29 sec 


( 1 ) Based upon uniform change over 20 minutes 


3.2.1 Star Tracker Errors 

The star tracker errors have been derived from analysis and test data 
associated with the single-gimbal engineering model tracker [2], The 
tracker and electronics errors are sunmarized in Table 3-3. 

3.2.2 Gyro Errors 

The development of gyro errors is based upon a gyro model having the 
characteristics as shown in Figure 3-2. The input axis rate, u, is gyro- 
scopically converted to a torque acting upon the viscously damped float. 
Angular motion of the latter about the gyro output axis is converted by the 
signal generator to an electrical output whose average value is essentially 
linear with respect to float motion. This output might simply be a contin- 
uous analog signal; or else a train of frequency modulated pulses, with the 
pulse rate dependent upon the float motion. The gyro electrical output is 
fed back to a magnetic torquer in order to rebalance the float. As indi- 
cated in Figure 3-2, appropriate selection of the gyro parameters usually 
results in either a rate gyro or a rate-integrating gyro. The former are 
generally characterized by fairly high bandwidths (e.g., 5 Hz or so) 
whereas the latter employ low bandwidths (e.g., 0.004 Hz) in order to effect 
their integration process. Clearly, in order to obtain angular outputs with 
a rate gyro, an additional integration must be performed. With frequency 
modulated pulse signal generators, this integration is easily accomplished 
for PAOS by use of a counter whose contents are read and reset to zero each 
integration period, e.g., 200 ms. Thus the gyro output used for PADS is a 
measure of the incremental change in angle measured about the gyro input 
axis over the (short-term) integration interval. 

Three basic uncorrelated gyro error sources have been postulated. 
Electronic noise can exist at the signal generator input (or equivalently, 
its output). Additionally, the various gyro mechanical imperfections give 
rise to a torque noise, n^H, acting directly upon the gyro float. Both 
these noise sources will be assumed zero mean and white. The fact that some 
self-correlation does exist in these noise sources is of little consequence 
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Table 3-3. Single-Gimbal Star Tracker and Electronics Error Characteristics 


SENSOR 


Electronic (Bias) 

.14 

arc-sec 

Thermal -Mech (Random) (30°C Range) 

.46 

arc-sec 

Noise Equivalent Angle 

,2 

arc-sec 

6IMBAL AND ENCODER 

Gimbal Alignment (Random) 

Runout 

.5 

arc-sec 

Perpendicularity 

.5 

arc-sec 

Gimbal Thermal -Mech (Random) 

Stress Relaxation 

.2 

arc-sec 

Thermal 

.2 

arc-sec 

Inductosyn (Random) 

6 

.7 

arc-sec 

256 0 

.5 

arc-sec 

READOUT ELECTRONICS 
Bias 

Zero Crossing Offset 

.13 

arc-sec 

Phase Stability 

.4 

arc-sec 

Random 

256 9 Excitation Harmonics 

.27 

arc-sec 

512 9 

- Cross Coupling 

.27 

arc-sec 

Gain Unbalance 

.13 

arc-sec 

A-0 Converter 

.13 

arc-sec 

Noise-Quantization 

0.3 

a re- sec 


] 


3-7 



FLOAT 

DYNAMICS 


SIGNAL 

GENERATOR 


DEFINITIONS 

= GYRO ROTOR MOMENTUM 

= GYRO DAMPING CONSTANT F 

= SIGNAL GENERATOR GAIN 

= TORQUE REBALANCE GAIN ^ 

= GYRO INPUT RATE 

= FIXED GYRO DRIFT RATE 

= GYRO OUTPUT 

= INTEGRAL OF GYRO OUTPUT FOR RATE GYROS 
= Cj/KK^ 

= WHITE SIGNAL GENERATOR NOISE 
= WHITE FLOAT TORQUE NOISE/H 
= WHITE TORQUE DERIVATIVE NOISE/H 


TORQUE 

REBALANCE 


RATE GYRO; K. 


H, xg « 1 


RATE INTEGRATING GYRO: K = C^/H, x » 1 


PULSE COUNTER USED IN 
RATE GYRO APPLICATION 


E[n^. (t)] = 0 


1 ,j = e,v,u 


E[n.(t^)nj(t2)] = 


o; 3 f 1 


Figure 3-2. Block Diagram of Gyro Model and Definition of Parameters 







as their bandwidths normally extend far beyond that of the gyro loop itself. 
Finally, there exists a gyro drift rate bias, by far the most significant 
error source unless appropriately compensated. While it is tempting to 
suppose that this drift rate holds constant at some value, b, it must be 
recognized that this bias may, itself, vary over the long time intervals 
normally encountered in attitude determination applications due to thermal 
and other effects. The model treats the drift rate as a random walk pro- 
cess, biased to the value b at some initial time, and wandering thereafter 
as a result of the zero mean white noise influence, n^. The fact that this, 
drift rate model implies an unbounded gyro output variance with time is of 
little concern here, as In^jl is numerically so small that for the typical 
time intervals indicative of attitude determination applications, a real- 
istic drift behavior is indicated. Experimental justification of this en- 
tire gyro model is demonstrated at a later point. 

The characteristics of the integral of rate gyvo output is summarized 
in Table 3-4. Noting the definitions introduced in Figure 3-2, Table 3-4 
shows how the integrated rate output, e, depends upon the various error 
sources previously discussed. Transfer functions are first presented re- 
lating these outputs to the error sources; and then the output standard 
deviation of each error is shown. Letting "E" denote expectation and g(t), 
the pertinent transfer function impulse response, the standard, deviations 
were computed by noting that 


ae'(t) 


[e^ = E { 

" t 


"t 1 

y'fi.j('r)g (t-x)dT 


f n|(x)g (t-x)dx ! 


i) 


3 ■ _ ' 



' t 

^•(A)j g(t-x)dx 

0 


dx 


= f g^(t-t)dT 

0 


(1 = e, V, u) 
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Table 3-4. Characteristics of the Integral of Rate Gyro Outputs 


GYRO ERROR SOURCE 

OUTPUT TRANSFER FUNCTION 
TO 

NOISE SOURCE 

OUTPUT STANDARD MVIATION AT TIME t 

— 

X > 0 

S 

LIMIT AS T V 

S 

(IDEAL RATE GYRO) 

WHITE FLOAT TORQUE 
NOISE/H 

1 r 1 



» 

DRIFT RATE DBRIVATrvE 
FLOAT TORQUE 'NOISE/H ' 

-irf ^ 1 


a tV2 
u 

• [\ H 

/2 

WHITE SIGNAL 
GENERATOR NOISE 

1 

T 

CTj s+1) 


D 

ONGERTAINTY- IN TEE 
DRIFT RATE BIAS AT 

t-0 ' 1 

1 

s(Tg S+1) 






is the standard dcrviaticm of the gyro bias uncertainty at t*0 























Note that only the signal generator noise in the rate gyro integral applica- 
tion results in a stationary process. This is because the noise is high 
pass filtered by the gyro loop dynamics, and thus the gyro output has zero 
noise power density at d.c. 


Assuming the rate gyro time constant to be small but f-nite, the sepa- 
rate effects of Table 3-4 can be root-sum-squared to yeild the dependence 
of the output standard deviation, to all the error sources: 


T . 0^ 

g e 
-Y- 


^ 2. . 2.2,1 2.3 

+ 0y t + 0^ t + y 0y t 


1 

2 


2 0 12 3 

It is seen that 0 contains terms dependent upon t , t , t , and t . It 

will be later observed that a minimum limit in 0 j^, dependent upon and 0 ^, 

exists if the drift rate bias is to be calibrated using the Kalman filter 

estimation. 


Figure 3-3 illustrates the results of data taken on two gas bearing 
gyros typically used in attitude determination applications - the Nortronics 
K7G and the Honeywell GG334. Fixing a specific sampling (or integration) 
time, t, the standard deviation 0 ^, of a large number of samples (with the 
drift rate bias removed) was experimentally established. This was done for 
integration times ranging from 0.1 second to 10 seconds. It can be seen 
that Figure 3-3 supports the model which was derived. For 0.1 < t ^ 10 sec- 
onds, the data indicates no strong time dependency and thus can be charac- 
terized as signal generator novse. For 10 < t ^ 2000 seconds, the slope of 
the data on this log-log plot is about 1/2, indicative of the noise source 
n^. As the sampling interval increases, the slope grows more steep, until 
finally it reaches 3/2, matching that of the drift rate variation noise 
source, n^. The analytical model derived for the data of Figure 3-3 is ex- 
pected to improve with use of current gyros. For example, based upon lim- 
ited data, 0 y = 0.015 appears a good choice for the Bendix 64-PM-RIG gyros. 
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Performance tradeoff curves were derived which relate the modelled 
gyro performance to the (attitude and gyro bias) estimation process. Error 
is normalized to (attitude update) sensor random noise. Figure 3-4 is a 
series of curves for estimated bias including gyro performance representative 
of the state of the art. Typically, 0.1 1 . _< 1.0 and 0.05 1 £ 0,5 for 

1 s^c £ ^ 5 sec and update period of 10-20 minutes. Estimates of gyro 

bias are on the order of 0.002 degree/hour - limited by the (ideal) case of 
a perfect attitude sensor: 


% = 




1 

T 


3.3 Covariance Analysis 

A PADS covariance analysis digital computer program was developed to 
evaluate PADS performance using the single-gimbal star tracker. The covar- 
iance analysis establishes a lower bound on performance under the essential 
assumption that the statistical models employed in the Kalman filter accu- 
rately represent the real world. 


3.3.1 Equation Development 

The equations are developed to establish the covariance associated with 

estimation of attitude and gyro bias. Thus the state vector, x, comprises 

three attitude variables (e.g., e.) and three gyro biases. The state error 
T ' 

covariance, E[xx ], is propagated according to the Kalman filter equations: 


^ j+1 " '^j+1 '"j ^j+1 ^j+1 


"j+1 


'"j+1 


Sj.,1 S .,1 . R.^^j 
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Performance Estimate Tradeoff Curves 






where: 


T 

P-- = E[xx ] at time t. just after an update 

0 w 

T 

S- = Etxx ] at time t. just before an update 

J J 

It is necessary to derive the state transition matrix, the state noise 
covariance matrix, Q, the measurement matrix, M, and the measurement noise 
covariance matrix, R, based upon the system model. 


The spacecraft is assumed to be three-axis stabilized to an earth 
pointing, local vertical reference frame rotating at orbit rate, The 
attitude of the spacecraft [b-framel relative to the reference [r-frame] is 
given by: 



PADS is assumed configured with gyro measured rates (attitude) in three or- 
thogonal axes [b-framel . The single-gimbal star tracker is mounted to pro- 
vide mechanical scan in a plane normal to the orbit plane (shown previously 
in Figure 3-1), with tracker axes [j-frame] related to the body axes [b- 
frame] as: 




~1 

0 

o' 


"b 


- 

0 

0 

-1 


^b 



0 

1 

0 


% 
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where 


X. = axis of gimbal rotational freedom 

y. = nominal sensor boresight at zero gimbal rotation 

J 

z. = form right-handed orthogonal frame with x. , y. 

\J J J 

From Section 3,2.2, each gyro is assumed to satisfy a mathematical 
model given by 


1 111 
03 ' = 03 + W + V 

m 


i = 1»2,3 


where is a white noise component; w^ is the random walk gyro drift rate; 
and 


03 = 6- - 03„ 0„ 

1 0 3 

2 

03 = 0o “ 03„ 

2 0 

3 

03 = 0- + 03 0- 

3 0 1 


Therefore; 


1 * 11 

0 = 0, - 03„ 0, + W + V 

m 1 0 3 

2 ' 2 2 
l3 = 0„ - 03 + W + V 

m 2 0 

3 * 3 3 

= 0-3 + “n 01 + W'" + V 

m 3 0 1 


For estimates: 




m 

2 

m 

3 

m 


01 - 0 ). 9o + w 
1 0 3 


02 - 


0o + 03. 0, + W"' 

3 0 1 
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Setting 60^ = 0^. - 0^ and 6W = w - w 


11 

602 - i ») Q 602 “ 6 w “V 

2 P 

60 g = - 6 W - V 

60 ^ “ -(»)q 602 “ <SW^ - 


1 . 


Also, as w is a random walk 


, 




= 0 


and therefore: 


• i * i • i i 

6w = w - w = u 


where u is random white noise. 

Since pitch and roll /yaw decouple kinematically, separate solutions can 
be obtained for pitch and roll /yaw. For pitch 


CD 

r 

2 

II 

[6W 

0+1 L! 


d+1 

1 


60 


+ 


/ 






2 


ds 


For roll /yaw: 


60 , 


*1 

6W 


60 . 


6W^ 


’ 0 

-1 

bi 

0 

o' 


60 


■-vl' 

0 

0 

0 

0 


6W^ 



-“o 

0 

0 

-1 


603 


-v^ 

0 

0 

0 

q 


6W^_ 


. 


6w''(t) = 6w''(t,) + u^' ds; i = 1,3; t,. < t < t 


j+1 
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Then for t ^ ^ t 



for which the honogeneous solution is developed as; 



t 

”cos W 

0 

Vl 

■■"j+l 

sin m T .,1 

0 J+1 

0 


■56j(tj)- 


_ 

0 


1 

0 

0 


sw^(tj) 



-sin u 

0 

'"j+1 

0 

cos T 

""j+l 





0 


0 

0 

1 




The (full) state transition matrix is summarized: 


$ = 


C 

0 

0 

0 

-s 

0 


-T 


j+1 

1 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 


-T 


0 

0 

d+1 

1 

0 

0 


where; C = cos S = sin 


S 

0 

0 

0 

C 

0 


0 

0 

0 

0 

"'^j+1 

1 
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The state noise covariance matrix for roll /yaw is developed from the 
particular solutions: 




B(t)/ 

0 


.t 


COSO) S 
0 

SCOSto s 
0 

-sino) s 
0 

-ssinu) s 
0 


r-v^ 



0 

1 

0 

0 




q' 

sinters 

ssinoj^i 

COSw s 
0 

ssinw s 
0 


-v^ 

ds = 

q' 

. 0 

0 

0 

1 






where B(t) is the matrix for the homogeneous solution. 
After very considerable algebraic manipulation, 


= E 


^ (1-COS 0)q t) 


+ T 


i[ 4 ] - E[qg = 


2 

= CT T 

u 




= o 


1-cos 0 ) T o 
0 Z 

_ ^ 


= e[Vi] = Ejq^qj] = £[,3,^] = j 

' E[q3qjl = E[q2q^] = E[q^q2] = 0 


Thus, the full state noise covariance matrix, Q, has the form: 


3-19 


3-20 




To establish the measurement equation, results previously developed 
for the two-gimbal tracker are used directly with the constraint that the 
inner gimbal is "locked". Recalling the measurement equation as: 



where 

yj is the sine of the inner gimbal 

y^ is the sine of the outer gimbal 

M M 

,a 2 are the star sensor outputs 

n^,n2 is the sensor noise 


and 



where u^ is the star unit vector defined in the reference axes as: 

The key constraint is applied by assuring that a^^ is first order , ag < 0 
and s a^.^ 1 
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Therefore: 


^1 "" ^1 ^2 ®3 " ^3 ®2 order) 


^2 ~ "®1 ®2 ^2 ^1 " ®3 


"3 = -®1 ®3 ■*• "2 ■" "3 ®1 


Substitution for u^. in the prior equation yields: 




- a. 


M 


+ n. 


+ d3 aj 


and 


y 



0 

^3 

■^2 




F;, M 1 

"®1 ■ “2 

= 





CD > 

ro 

+ 



CO 

_51 

0 

0 

- 


h. 


a + a M 

_ ®2 ®3 “l _ 


5y 


■ 

0 

^3 

-^2 


"2 

y-y = 

"3 

0 

0 

66 + 

:^3 '^1. 


Therefore, the measurement matrix, M, is given as 
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The measurement noise covariance matrix, R, is given as; 




3.3.2 Performance Eval uati on 

The covariance analysis digital computer program was used to evaluate 
performance potential and tradeoffs. Those elements evaluated for tradeoffs 
include star sensor noise and gyro random drift characteristics. The com- 
puter program provides a time history of the error covariance. Summary data, 
used as the performance measure, includes the average estimation error covar- 
iance and final value of gyro bias estimation error covariance. 


The program was run with several sets of initial conditions to evalu- 
ate both the convergence properties and steady-state performance. Conver- 
gence to steady-state accuracies from initial attitude uncertainty of 
0.1 degree (I 0 ) and gyro bias uncertainty of 0.1 degree/hour (I 0 ) was 
achieved within four (4) star sightings. Cases for evaluation of steady- 
state performance had initial conditions of 10 arc-seconds (la) in attitude 
uncertainty and 0.05 degree/hour (la) in gyro bias uncertainty. Small ini- 
tial transients were noted in the attitude covariance. 


Table 3-5 summarizes the observed covariance analysis results as sum- 
marized by the computer program performance measures. The gyro parameters 
are consistent with a range of values based upon observed test data (see 
Figure 3-3), The star tracker errors are consistent with those developed 
in the error analysis of Table 3-3. Evaluation of performance goal of 0.001 
degree (= 3.6 arc-seconds) is achievable in the high altitude (geosynchronous) 
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Table 3-5. Table of Performance Results^^^ 


Trackei^oise 

CTp (^) 

Gyro Random Drift 
(sec/sec) (sec/sec^) 

(rad/sec)^^^ ^ (rad/sec)^^^ 

Attitude ^imate 
a (sec) 

®i 

Gyro Bias Estimate 
CT. (deg/hr) 
i 

2 

.115 

3.33x10"^ 

4.6/6. 1/6.2 

.0024/. 0026/. 0025 

2 

.05 

3.33x10“^ 

3. 2/4. 2/4.1 

.002/. 0022/. 0021 

2 

.05 

1.6x10"^ 

2. 6/3. 1/3.1 

.0011/. 0012/. 0012 

3 

.115 

3.33x10“^ 

5.0/6.4/6.5 

.0024/. 0026/. 0025 

3 

.05 

3.2x10'^ 

3. 5/4. 4/4. 3 

.0016/. 002/. 0018 

3 

.05 

1.6x10"^ 

3.0/3. 5/3.4 

.001/. 0012/. 0011 

3 

.015 

3.2x10“^ 

3. 1/3. 8/3. 7 

.0014/. 0018/. 0016 

3 

.015 

1.6x10'^ 

2. 5/2. 7/2. 7 

.0009/. 001/. 0009 


cases - geosynchronous orbit, real stars brighter than m^ = +3.5, 
star tracker boresight South, gimbal freedom ±45 degrees 


application using the current PADS single-gimbal star tracker and available 
state-of-the art strapdown gyros. Certain general conclusions may also be 
developed, namely: 

f Small star tracker noise variations are observed to have little 
effect on overall performance 

• Gyro drift characteristics have the largest effect on performance 

• Changes in the gyro drift term, have a moderate effect on 
attitude and bias estimation 

• Changes in gyro drift term, have (almost) a direct effect on 
the bias estimation 

A typical time history of the attitude estimation performance is shown 
in Figure 3-5. The time-history plotted is the square-root of the trace of 
the attitude error covariance matrix. Selected case history output data is 
included in Appendix B. 
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appe;ndix a 

COMPUTER PROGRAM FOR STAR AVAILABILITY 


A-1 


A.l Program L00K1 

Program L00K1 is a FORTRAN program to determine and order stars ob- 
served by a star sensor mounted on an earth orbiting spacecraft. The 
program works from stars in a general catalog (input file TAPEl) which in- 
cludes RA, DEC, magnitude, and catalog number. Each star is tested in 
turn for the position in orbit that observation may be made including tests 
for sun interference, earth interference, and constraints of orbit, time of 
year, and sensor/spacecraft geometry. If observation is determined, the 
star is ordered in a v/orking catalog on the basis of orbit angle from as- 
cending node at the position of observation. The star LOS is defined by 
orbit angle and angle measured from zenith at the indicated orbit angle. 

The input file used for current investigations is given in A. 2 and the 
output file of ordered stars used for covariance analysis given in A. 3. 


Jinioa ^PRDGRAM LOOKiCINPUT .OUTPUT, TAPFi=INPUT5 TAPE?.j.TAP£3| 

00110 C TAPE3 CONTAINS MASTER STAR CATALOG 

00120 DIMENSION NNS { A 00 1 , NUS (40 0 ) , GAM (4 0 0 ) , ALPH ( 400 ) 

_ REAL INCL,.MV,MVtIM 

|0014Q REWIND 3 

■00150 DATA OMEG.INCL, S/0.,0.,0./ 

J]iLL&D DATA HVLIM».SUNL IM-.BETA, GAMLTM/3.5.45. .-45. . 45./ 

00170 NAMELIST/DATA/OMEG,INCL,S,MVLIM,SUNLIM,eETA,GAMLIM 

00180 READdtDATA) 


DD19.a DIS PL A Y ♦ OMP ,.I NCL-j. S ~ JLi-OME G 

00195 display »HVLIM,SUNLIM,BETA,GAMLIM= MVLlM,SUNLiM,BfTA,GAMLIM 

00200 OEGRAD=3. 1415927/180. 

IQQ.21 n ._-,RADDEi;=la/DEGRAIl_.. 

00220 SUNLIM=DEGRAD*SUNLIH 

00230 GAMLIM=DEGRAD*GAMLIM 

,011. 24 0 _QME.G=DEGRA.D.T11MJELG. 

■00250 INCL=DEGRAD»INCL 

IDD260 S=DEGRAD*S 

itm.2.7.Q_ _ACCEPTJ3)...ilS_ 


00280 3> L=1 

00290 do DO 500 1=1, NS 

aai,ILQ ACCEPT O) ..NSL,.NSC,PA,DEn»m.. 

00310 C TEST ON STAR MAGNITUDE 

00320 IF<MV.GT.NVLTM)G0 TO 500 

‘ 00330 C FORM STA„R_.G.IRECTI0N COSINES . 

00340 RA=DEGRAD*RA 

00350 DEG=OFGRAD*DFC 

00360 A1 = C 0 S C DEC ) W S(RA 5 

100370 A2=COS!DEC)*sfN(RA) 

'00380 A3=SINCDEC) 

iim.33.Q C, TEST., SUN .INTERFER£J1Q£, 


00400 

00410 

MAZE 


Sl=C-3SfS) 

S2=STN(S)»COSC.408) 
JS3=SIN IS) *SIHLa-A]lSJ 


100430 D0TSUN=A1»S1+A2»S2 + A3*S3 

00440 IF(DOTSUN.GT.COSCSUNLIM))GO TO 

ittOAEQ NNS IU^ 


500 


00460 FLAG=D. 

00470 C DETERMINE REF PLANF-ORBIT PLANE INTERSPCTTON 

00480 _ Bin SXN (I NCLTYS I N (OME G) 

|00 490 B2=-SIN ( INCL) ♦COS ( OMEG ) 

jOOSOQ B3=G0S{INCL) 


'Ss. : i 

; r 

!oQ510 

00520 ¥2=A3*B1-41*B3 

00530 V3=fll»B2-A2*Bl 

00 540 H1=R2»V3-B3»V2 

p550 W2=B3*V1-31»V3 

100560 H3=B1*V2-B2*V1 

’ 00570 if=SQBT (W,1*H1+W2.^H2+H31H3)_ 

00580 W1=W1/W 

00530 M2=W2/W 

00600 J3=H3/W_ . .... 

‘00610 C COMPUTE ANGLE FROM ZENITH FOR STAR LOS 

00620 GAMMA=AC0S(W1*A1-I-W2»A2+N3*A3) 

00.630 IFf (Al*Bl-t-A2*B2-»-A3»B3).GT.Q.)G0 T O 10 0 

00640 GO TO 200 

00650 iOO GAMMA=-GAMMA 

^660. ..2.0.0 _I.F C A_8S.CGAMMA-nEGPA.P.^^ET.A ) . GT_. GAM.UIM1LG.O... TQ. 5.0.0_. 

00670 GAM{L)=RAOOFG»GAMMA 

0068 0 C COMPUTE ORBIT ANGLE FOR REF PLANE INTEP^SECT 

00.650 Aill=C.OS(OM.EG) .... 

00700 3> AN2=SIN(0MFG) 

00710 SALPHA0=83» (-AN2»W14AN1»W2) + {31»AN2-B2’^AN1)*W3 

00720 _ CALPHA0=AN1»W1 + AN2*H2 _ 

'00730 ALPHA 0=AT AN2(SALPH AC, CALPHAO) 

00740 ALPH(L»=RADOEG*ALPHAO 

;Q0.7.5.0 . . . If.CA_LPH(L)..LT,0.),ALPH.(L)=.ALJ^.H(L.)+36 

6o760 C SORT ON BASIS OF ORBIT ANGLE 

00770 250 NUS(L)=L 

M 78 0 I TFMP 2=.N_U _S { 1 ) 

00790 00* 300 J=i,L 

00800 M=NUS(J) 

io OeiD IF(ALPH(L) ■LT.ALP HCM ))GO TO 400 ; 

00820 ITEMP2=NUS(J+1) 

00830 300 CONTINUE 

|ULaM.__„4.aO .NUS (Ul=JiUSJLU. 

100850 DO 450 K=J,L 

j00860 ITEMP1=ITEMP2 

00 870 IXE MP.2= N U S ,( K+1 ) 

00880 NUS CK+ii-ITEMPi 

00890 450 CONTINUE 

0090 0 c comple.tp sort 

b0905 C TEST FOR ALPHA+ISO 

B0910 IFtFLAG.GT.l.lGO TO 480 


01020 
'01030 C 
iD lQtf fl, C.. 
01050 
01055 

QIO&Q 

|0107 0 
101080 

01100 


FlAe=2 .., . 

IF(ABS (180.-GAH(L)-BETft) .GT.GAMLIM* 
L=L+1 

...NMS LL).=J: - 

GAHtL)=18Q.-GAMfL-l) 
ALPH(L)=ALPH(L-1) +180. 

GO TO 250 
CONTINUE 

._L=„L+i. ^ 

CONTINUE 

END OF SOPT ITERATION LOOP 

„?SIEJ_ilROER£a-.G-flIALJlG. 

L=L*1 

DISPLAY (2) L 

. JIO JJin. , . 1 = 1,.L 

J=NUS(I) 

DISPLAY <2) T,NNSCJ) ,ALPH(J) ,GAM(J) 

.CONTINUE 

END 


GT.GAMLIM»RADDEG)Gn TO 



A. 2 Yale 
Data 

• 

t 


Catalog of Bright Stars to 4th Magnitude 
includes (from L to R): 

Working catalog order (ordered by magnitude) 
Ya^'e Catlog Number 
Right Ascension 
Declination 




Magnitude. (S-20) 


2491 

100,91355 

-16.67130 

4210 

160.93650 

-59.51870 

._2.32j6u 

9*5 .7 991*5 

-52.66630 

7001 

273.94600 

38.74930 

1713 

78.22530 

-8.23967 

AZ2_ 

24.10905 

-57.42000 

5267 

210.35700 

-60.20230 

1708 

78.54330 

45.96600 

■■2943 ■ 

-_lJL4,.5Ai}lO 

5.31833 

5056 

200.84850 

-10.97430 

534 0 

213.52800 

19.35900 


_L2. -ZSfLk.. 

13 5i»59 

14 4853 


I 16 3982 

I 17 4730 

19^ 2618 
20^ 2061 
21 7924 

' 22 6527 

23 1790 

i- 2h iaJL3_ 


2.9.7^. 3JL2DQ. 

219.32250 
191.43600 
_ 297^6I15.I1_„. 
151.63950 
186.17700 
.. 343.-9Jt4JUl. ^ 
104.32395 
88. 33125 
.-3UL..J16Ma._. 
262.82400 
8 0.8 2720 

_ _a.3.6153.,Q „. 


32. 83JL2JL 
■60.69170 
•59.51300 


12.13100 

-62.91300 

-2a.jrS8JLiL 

-28.92130 

7.39433 

45«,1A72J1. 

-37.07730 

6.32167 

-JL*22.257_ 


31 5191 

32 6134 

^..3? 3£S5- 

I 34 7790 

! 35 2294 

49J15 . 

37 2990 

38 6879 

_ 39 4 731 

40 2004 

' 41 7121 


„ 1JL3 aJ735,. 
206.54850 
246.83250 
.X3JL...2JL655_ 
30 5,73750 
95.30100 
_193,.l_34il0. 
115.80 930 
275.478 0 0 
..i8.6^1J_150__ 
86.53785 
283. 290 00 


49.48670 

•26.35970 

•69 .57500 

•56.84100 

■17.93300 

56.1 3700 

28.10170 

•34.40030 

•6 2.91300 

-9.67800 

■26.33970 


-1.4600 
- 1. 0000 
-.5754 
. 0400 
. 0487 
.2599 
.3440 
.6530 
.6880 
.7019 
.7217 

.8356 
.9698 
- .-97-7-2_ 
1.2420 
1.2480 
■ l.,2490 
1.2540 
1.3320 
1,-3.490 . 
1. 3500 
1.3820 


25 

1457 

68.49180 

16.43200 

1. 5000 

26 

1948 

84.76245 

-1.96700 

1. 5060 

27 

-1231... 

81. 03390 

28.57170 

1. 5090 

28 

3207 

122.12130 

-47.24300 

1. 5250 

29 

8425 

331.52250 

-47.13100 

1.5890 


1.6110 
1. 5380 
1.6380 

.. It ZQ.l O 

1.7080 
1.7100 
1.7390 
1.7300 
1.8090 
1. 8320 
1.8 420 
1.8780 



kfu ; ■ 

'f 

» 

i 42 15.,_ 1,65666 ..2.8,89630 lj..913L 

43 2421 98.93760 16.42830 1.9300 

44 2088 89.26005 44.94430 1.9310 

45 316 5 120.59685 -39 .90370 1 .95 50 

[“46 1852 ' 82.56930’” -.32267 1.9660 

47 3485 130.94400 -54.59200 1.9910 

48 59J58 239.52000 26.01300 2. 0000 

’"49/ 2891 113.10735 ’31.95700“ 2.0210 

50 5132 204.43200 -53.29100 2.0420 

\S_1 5469 219.91800 -47,25270. 2. 05 40 

^ 52 936 46.49190 ” 40.81970 2.0370 

‘ 53 5440 218.33550 -42.00270 2.1280 

... 54.. .481.9 189. 9 1 0_5 0 - 48, 77970 2. 1390 

55 6580 265.03350 -39.01630 2.1520 

56 1017 50.47560 49.73100 2.1920 

57__ A953. 239.58000 -22.520 30 2.2020 
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211.11150 

-26.51900 

3.9020 

-3.57 - 

5603- 

- 22 5. 5-19 0 Q_ 

... -_2 5» 147,3.0. 

... 1 » 9Qsa . 

358 

5993 

241.20300 

-20.57600 

3.9080 

359 

5747 

231.60450 

29.21900 

3. 9090 

u3.6Q. . 

6324 

254.74950 

.-3,0_,.g677.0„ 

_ 3,910 0. 

361 

4802 

188.96100 

-48.34630 

3.9100 

362 

7852 

307.89750 

11.18670 

3.9110 

-363 . 

--2.19 

11.76075 

_ 57.63530. 

3.9130 

364 

6623 

266.28300 

27.75 0^0 

3,9140 

365 

5744 

231.04350 

59.08570 

3.9180 

t.366- 

J652 

144.83416 . 

1,0_, 05300.. 

3.9190 

367 

7590 

299.16450 

-72.99070 

3.9190 

368 

2216 

93.20385 

22.51130 

3.9190 

369 

1251 

6 0. 335.56 „ 

5.90367 

3.9200 

370 

8597 

338.40150 

-.29233 

3.9240 

! 371 

7420 

92.21050 

51.65970 

3.9260 

L^72-_ 

-332-3. 

12. 6^8-6 115 

60.83000 . 

3. 9430 

373 

1336 

63.49515 

-62.55170 

3,9450 

374 

1949 

84.76245 

-1.96700 

3.9520 

3Z5._ 

.4-826.- 

-- 18 9.9855J... 

-1.263 00 

3.9530 

376 

3757 

142.2i320 

63.21400 

3.9540 

, 377 

8278 

324.55050 

-16.81970 

3.9560 


I 


I 




378 

1756 _ 

79^5006.1 . 

._r.13-2J.730_. 

3.9600 

379 

80 28 

313.97850 

41.03630 

3.9600 

380 

1922 

83.33055 

-62.5060G 

3.9680 

381 

580 

30.13155 

72.25230 

3.9700 _ 

382 

6698 

269.28600 

-9.77767 

3. 9700 

1 383 

2657 

109.55350 

-15.58230 

3.9710 

384 

7234 

28 6.20300 

-27.71770 

3.9780 

385 

6095 

245.10750 

19.22930 

3.9810 

386 

8465 

332.41800 

58.03000 

3.9840 

387 

.'•3_86_. 

169. 84500 

6.22033 

3.9860 

1 388 

269 

13.71575 

38.31300 

3.9880 

, 389 

1273 

61.54905 

47.62600 

3.9890 

■ 390 

_6714 _ . 

269.73300 

2.93333 

3.9900 ^ _ ^ - 

391 

921 

45.75045 

38.70300 

3.9950 

392 

6927 

275.68350 

72.71630 

3.9950 

393 

1983 

85.76250 

-22.46130 

3.9960 _ 

394 

20 85 

8 8.71330 

-14.17230 

3.9960 

395 

5367 

214.62000 

-37.74130 

3.998 0 

^ 396 

5987 

_241 ?M.750 „ 

-36.70930 

4.0000 





V.; 3r 

A. 3 Ordered Working Star Catalog 

Shown is the ordered working catalog derived from Program L00K1 for 
single-gimbal tracker, boresight South, ±45° FOV gimbal , m^ = +3.5. Data 
output includes: 

• Local catalog order (ordered by. orbit angle measured from 
ascending node) 

0 Working catalog number (see A. 2) 

0 Orbit angle of observation 
' 0 Angle from zenith 

Note: Output is provided on file named TAPE2 


A-17 


:.1M _ 

1 «7 1.647005318 1.2947D3E+02 

2 90 3.336005318 1.21437E+02 

3.. .. 182 ^ 5. 99126; „7. 744?ZE+ 01 

4 17 6.177005318 1.17087E+02 

5 39 6.181505318 1.17087E+02 

L 6„. , „7 , 320 0 0 531 8__. 1.23076E >02 

7 79 8.788505318 1.11Q537E+02 

8 54 9.910505316 1.312203E+02 

9 120 1.10i(9J0 5_32E+Gl. 1.1207035^02 

[To 14 1.143&6 o532E+01 1.20487E + 02 

! 11 6 2.410g05E + 01 5.742E+01 

’ 12 50 2. 443200532E 4^01 1 .26709E>02 

13 64 L 83'5300532E + 01 1.3287E + C2 

14 157 2.942385E+01 6.1731E+01 

^15 _ 7 3.0 35700532E-»- 01 1.197977E+02 

[lE 213 3.3821855+01 5.16753E+01 

; 17 195 3.430650532E+01 1.34Q92E+02 

_18 ^13 3.93225D532 E+01 1. 1930835+02 

193= 66" "3.‘g322505'32E+0l l.’"ig3083E+02 

20*^ 51 3.9918005325+01 1.3274735+02 

_21..“ . 206. .„.3.9«375.0532E+J.l, 1.15 18075 + 02 

22 127 4.8934505325+01 1.1144135+02 

i 23 156 5.8038005325+01 1,1667435+02 

165 6.A31 5855+01 _ 5 .51 185+01 

25 92 >.’1265 005325+01 1.110295+02 

26 231 8.0617505325+01 1.2450075+02 

27 .169 8.063 4005325+01 1. 2365075+02 

28 110 8.2305'D0532E+Oi 1.301455+02 

i 29 3 9.5799155 + 01 5. 266635 + 01 

1^0.. 1J8 9.6115505325 + 01 1.3401635 + 02 


31 

220 

1.0196225+02 

6.189935+01 

32 

180 

1.1897915+02 

5.289275+01 

33 

28.... 

1.2212135+02 

4.72485+01 

i 34 

84 

1.2545495+02 

5.939235+01 

f 35 

34 

1.2573750535+02 

1.231595+02 

l36 

208 

1.2983135+02 

5.279775+01 

37" 

137 

1.29876455+02 

4.65315+01 

38 

47 

1.309445+02 5 

.45925+01 

39 

177 

1.37517755+02 

5.88255+01 

i 40 

33 

1.38206555+02 

6.95755+01 

! 41 

72 

1.39046855+02 

5.91255+01 




Cf >~ lo) 





i ■ • 

V 


'-42.- 

63- - 

-.±.-4.0 265 6E*il2.-- 

5L.4J3,693E+iLl 

43 

207 

l,4656245E+02 

6.4908E+01 

44 

217 

1.4891775E^-02 

5.4419E+01 

,-45_ 

. 29 

■ 1.51 5??5053F-<-n? 1.2266QE+J12 

I 46 

178 

1.532355E+02 

6.98633E+01 

' 47 

205 

1.577025E+02 

6,lF077E-»-01 


M im . . _ JU6.CL15 as a5.3,E+iL2 UJi23.klE±JiZ 


49 85 1.60434E+02 6.42077E+D1 

50 2 1.609365E*02 5.95187E+01 

^3.1 JL85 - l^,&13265E+02 4..a£52£±01___ 

: 52 152 1.73553E+02 6.28297E+01 

53 87 1.316470053E+02 5.05297E+01 

.54. an .. 1.833360053E+02 5.85 63Ef01 

55 182 1.8599126E-I-02 1.025573E+02 

56 17 1.861770053E+02 6.2913E4-Q1 

_57 3.9„„ l.A61815a53.E+J}.2 _ 6...,29 jL3.E+JL1 

‘58 62 1.373200053E+02 5.6924E+01 

f 59 79 1.887885053E+02 6.89463E+01 

IfiJEL 54 - 1.89ai.05053£+02 .. A,37797.E*D.t.. 

61^ 120 1.91049Q053E+02 6.79297E + 01 

62 14 l.gi4360fl53E+02 5.9513E+01 


-63- , 

...6 

2.Q4109D3Eta2 .1 

.22SflE*.l).2 

64 

50 

2.Q44320053E+02 

5.3291F+01 

. 65 

64 

2.083530053E+D2 

4. 713E+01 

;-6j6_ .. 

157 

2..0 9.423a5.EtQ2 . 

l,18269£+02 

67 

7 

2.103570053E1-02 

6.02023E+01 

68 

213 

2.1382185E+02 

1.283247E+02 

-69 

. 195 

^ -2..14356.5fl_5-3FtA2 

4.5908E+01 


70 13 2.193225053E*02 6.06917E+01 

i 71 66 2.193225053E+02 6.06917F401 

LZ2.. ..5JL_. ..E^iagxa 0053E-t:0 2 4.72527E-»-01 


73 206 2.199375053E+02 6.48193E+01 

74 127 2.289345053E+02 6.85587F+Q1 

-25. 15 5 __ 2.i^38J3J.8 Q Q53E^-02 6.33257F4Q1 

: 76 165 2.483t585E + 02 1.24882E+Q2 

77 92 2.512650053E+02 6.8971E+01 

\ 78 221_ . 2.606175053 F*0 2 5.54993E-t-01 

79 169 2.606340053E+02 5.63493E+01 

80 110 2.623050053E+02 4.9855F+01 

81 .3_.._2._7579915.E+0 2 1.2 7333 7E »02 _ 

•82 188 2.761155053E+02 4.59837E+01 

, 83 220 2.819622E+02 1.181007E4-02 



t- 


I ; 

! 

LM_,.J.9a -2..989791Et.02 l,.271Q73E+0_2 .. 

35 23 3.021213E+a2 1.32752E+02 

86 84 3, 054549E1-02 1.206077E+02 

,.^7 3_4 3 _. 057375g53E+g 2 5.6841Et01__ 

; 88 208 3.098313E+02 l,272023E+02 

89 137 3.0987645E+02 1.33469E«-02 

• 90 47 . 3,10944E+02_ 1.254D3E ^-02 _ _ 

91 177 3.17S1775E+d2 i.21175F+02 

92 33 3.1820655E+02 1.10425E4-02 

93 72. 1. 19D468_5E+02 1. 20875E»02 

94 63 3.202656E + 02 ’ 1.251307E+G2 

95 207 3.2656245E+02 1.15092E+02 

_96 217 3 .2.891 775E+ 32 1.25581E»02 

97 29 3.3i5225d53E+02’ 4.7131F«-01 

98 178 3.332355E+02 1.101367E+02 

99 2 0.5 .3,.3.7.7 0 25E+ 0 2 1. 184923E»02 

100 106 3.4015950535+02 4.7059E+01 

: 101 85 3.40434E+02 1.157923E+02 

' 102 2. 3 titO 9 3 6 5E+.0 2 _ 1.. 20 k3 13E+02 

103 :t» 185 3.413265E+02 1.30748E+02 

104^152 3.53553E+02 1.171703E+02 


APPENDIX B 


COMPUTER PROGRAM FOR SINGLE-GIMBAL 
PADS COVARIANCE ANALYSIS 


B.l Program PADSC0 

A FORTRAN digital computer program, Program PADSC0, was developed 
to implement the covariance analysis algorithms derived in Section 3.3. 
The program operated using star data from file TAPE2 and input data de- 
fining orbit and filter characteristics. A source program listing is 
provided, followed by a summary of output data in B.2. 



I 0010 0 
00105 
00110 
00120 
00121 
00125 
00130 
00140 
00150 
00160 
00170 
0018 0 
00190 
00200 
00210 
00220 
00230 
0024 0 
00250 
00260 ' 
00270' 
00280 
00290 
00300 
00310 


00320 

00330 

0034 0 

00350 

0036 0 

0037 0 
00380 
00390 
00400 
00410 
00420 
00430 
00440 
00450 
00460 


00470 

00472 


PROGRAM PADSCOd NPUT, OUTPUT, TAPFl. TAPE?. TftPF5=:THPlJT> 

TAPE2 CONTAINS STAR DATA 
DIMENSION P(6t6) tS(6,6) ,SS(3 t2) 

DIMENSION ALPH(400) ,GAM(400»,N(400) , NNS (4 00 ) ♦ A ( 3) 

COMMON HO 
REMIND 2 

NAMELIST/PARAM/Tl.T2tT3,H «SN4SV«SU.H«TF«TG 
DATA T1,T2,T3/10.,10.,10./ 

DATA W/,05/ 

DATA SNtSV^SU/3.« ,115t»333E-4/ 

DATA H, TF,TG/19300, ,86400. ,2160 0. / ~ 

REA0(5,PARAM> 

WRITE (1,2) Tl,T2«T3tH,SN«SVtSU»H 
PRINT 2,T1,T2,T3,W,SN,SV,SU,H 

FORMAT (X,11HT1,T2,T3 = ,3 (F6. 1, 2X> /X, 4HW = ,E9. 2/X, IIHSN, SV, SU 
13(E9.2,2X) /X,4HH = ,^7.1/) 

WRITE (1,3) — — 

FORMAT (4X,1HT,7X,2HS1,6X,2HS2,6X,2HS3,6X,2HPS,6X,2HRS, 
16X,2HTl,6X,2HT?f 6X , 2HT 3 , 6X , 2HPT , 6X , 2HRT/) 

ACCENT (2) NSTR ~ 

ACCEPT (2) (N(I),NNS(I) ,ALPH(I) ,GAM(I) ,1=1, NSTR) 

DEGRA0=3. 14 15927/180. 

IR=0 “ ^ ^ 

INITIALIZE TIME 
T=T9=0. 


INITIALIZE PERFORMANCE MEASURE 
00 11 1=1,3 
SS(I,1)=SS(T,2)=0. 

INITIALIZE ERROR COVARIANCE MATRIX 
DO 10 1=1,6 
DO 10 J=l,6 
P(I,J)=0. 

P(1,1)=T1*»2 

P(3,3)=T2**2 

P(5,5)=T3»*2 

P(2,2)=P(4,4)=P(6,6)=M**2 
H0=SQRT (1.875E5/3./(3440. ♦H)**3) 

RT=SQRT (Tt*»2+T2*»2+T3»*2) ~ ~ 

RS=SQRT (Sl*»2+S2*+2+S3**2) 

PS=SDRT (Sl**2+S?**2) 


pT=SQRT(T1**2+T2*»2) ^ 

PRINT 8,T,RS,RT 



00760 C INCREMENT TIME TO NEXT STAR 

00770 30 T=T+OT 

00780 C TEST FOR END QR RUN,TF 

00790 IFCT.GT.TE) GO TO 100 

00800 CALL FILTER(P,S,nT,SN,SV,SU,A) 

00810 S1=SQRT(S(1,1) ) 

00820 S2=S3RT (S(3,3) ) 

00830 S3=SQRT(S(5,5) ) 



00840 


00850 
00860 
00870 
00880 C 
00890 
00900 
00910 
00920 
00930 
00940 
00950 
00960 
00970 
00980 C 
00990 


01000 
01010 
I 01020 


01030 , 
01040 
01050 
01060 G 
01070 C 
01080 C 
01090 
01100 
OHIO 
01111 
01120 C 
01130 
01140 
01150 
01160 
101170 C 
I 0118 0 
0119 0 


01200 
01210 C 
01220 


01230 

01240 


Tl=SdRT (P(l,m 
T2=SaRT(P(3,3) i 
T3=SQRT(P(5,5) ) 

GO TO 20 

COMPUTE EOR PERP. MEASUR£:MEAN, SIGMA 
100 CONTINUE 

DO 50 1 = 1^3 

S3 (I,1)=SS (T,l) /T9 
SS (I,2)=SQRT(SSCI»2)/T9-SS{I,1) **Z) 

50 continue 

PRINT 5 

5 FORMAT (5X,4HME AN, 4y»5HSIGMA/) 

DO 60 1=1,3 

60 PRINT 6,SSCI,1) ,SS(I,2) 

ESTIMATED GYRO BIASES 
PP3=SQRT(P (2,?) ) 


PP2=SQRT{P (4,4)) 
PP3=SQRT(P(6, 6) ) 
PRINT 7 ,PP1,PP?,PP3 


7 F0RMAT(X,14HGYR0 BIASES = , 3 { El 0, 3 , 2X ) ) 
6 F0RMAT(2(?X,P6.2) ) 

END 


SUBROUTINE FILTER ( P ,S , OT, SN , SV, SU , Y ) 

DIMENSION P(6,6),S(6,6) ,Q{6,6) ,R(2,2) ,RH(2,6) , PHI (6,6) ,A(6,6) , 
1B(6,2) ,C(2,2) ,0(2,2),G(6,2) ,H(6,6) ,Y(3) 

CCMMON WO 

ZERO MATRIX ELEMENTS 

DO 10 J-1,6 ^ 

RM(1, J)=RM(2,J)=0. 

DO 10 1=1,6 

Q(I, J)=PHI(I, J)=0 . 

MEASUREMENT NOISE COVARIANCE MATRIX 
R(1,2)=R(2,1) = 0. 

R(l,l)=(SN/2.)*»? 


R(2,2)= (SN*Y(3) )*^2 

STATE TRANSITION MATRIX, PHI 

DO 20 1=1,6 


PHT(I,I)=1. 

S8=SIMIWo^0T) 



01241 



C8=C0S(Hn»DT) 



01242 



PHI(1,1)=DHI(5,5)=C8 



01243 



PHT(i,5)=SS 



01244 



PH1(5,1)=-S8 



01245 



PHI(1,2}=PHI{3,4)=PHI (5,6)=-nT 



01250 

C 


STATE NOISE COVARIAhCE MATRIX, Q 



01260 



T4=4.»DT**3/3.-2.»DT» (1 .-C8) /W0»^2 



01261 



0(3,3)=DT*SV*»2+SU»^2»DT»+3/3. 



01262 



0(1,1) =0 (5,5) =DT»SV»»2+T4»SU*»2 



01263 



Q(1,2)=Q(2,1)=0(5,6)=Q(6,5)=( (1 . - C8 ) / WO »» ?-OT»*2 ) 


01264 



Q(1,6)=0(6,1)=Q (2,5) = Q(5,2)=SU-^*2»(0T-S8/W0)/W0 



01265 



0(2,2)=Q(4,4)=Q (6,6)=DT»SU»»2 



01266 



Q(3,4)=Q(4,3)=-(SU*CT)»»2/2. 



01280 

c 


MEASURFMENT MATRIX, RM 



01290 



RM(1,3)=RM(?,1)=Y(3) 



01300 



RM(1, 5)=-Y(2* 



01310 

01320 

c 


PROPAGATE ERROR COVARIANCE, S=PHI»P»PHIT+Q 
DO 30 1=1,6 

■'/ isfi 


01330 



DO 30 J=l,6 

. . .. 


01340 

CO 

1 


A(I,J)=0. 

2 


01350 

cn 


DO 30 K=l,6 

- ^ 


01360 


30 

A( I, J) =A (I, J) fP (I,K)»PHI( J,K) 



01370 



DO 40 1=1,6 

'.i 


01380 



DO 40 J=l,6 


01390 



S( I,J) = Q(I,J) 

.. 


01400 



DO 40 K=l,6 



01410 


40 

S(I,J)=S(I, J)+PHI(I,K)»A(K,J) 


01420 

c 


COMPUTE OPTIMAL GAIN,G 

;i'. hR ...... 


01430 



DO 50 1=1,6 

ISJ 


01440 



DO 50 J=l,2 


01450 



B(I,J)=0. 




01460 DO 50 K=l,6 

01470 50 8(I,J)=8(I, J) +S(I,K)*PH (J,K) 

01480 on 50 1 = 1,? 


01490 DO 60 J=l,? 

01500 C(l, J) 

01510 DO 60 K^l,6 

015?0 60 C(I,J)=C(I,J)+RM{I,K)»a(K,J) 

01530 DEL=C (1,1)*C<2,2)-C (1,?)*C (2,1) 


01540 D(1,1)=C(?,2)/DEL 

01550 0(1,?)=-C(l,?r /DEL 

01560 D( 2,1)--C(?,1) /DEL 


01570 0f?,2)=C(l,l)/DFL 

01530 DO 70 1=1,6 

01590 00 70 J=l,2 

01600 G(I,J) = 0. 

01610 90 70 K=l,2 

01620 70 GtI,J)=G(I,J)+B(I,K)»D(K,J) 

01630 C UPDATE ERROR CQ V APIANCE,P 

01640 DO 30 1=1,6 ' 

01650 DO 81 J=l,6 

01660 81 H(T,J)=-6(I,1)*PM(1,J)-G(T,2)’^RM(2,.J 

01670 80 H{I,I5 = H(1, 1)4-1. 

01680 00 90 1=1,6 

01690 00 90 J=l,6 

01700 A(I,J)=0. 

01710 DO 90 K=l,6 

01720 90 A(T,J)=A (I,J) 4-H(I.K)»G(K, J) 

01730 C ASSURE P03 . DEF. P 

01740 00 100 1=1,6 

0175 0 DO 100 J=l,e 

0176 0 oo 100 Pa,J»=CA(I,J» 4-A{J,i))*.5 

01770 RETURN 

01780 END 




B.2 Summary of Test Cases 


i 


Star Data Used From Appendix A 
Geosynchronous Orbit 

Tracker - Nominal Boresight S, Gimbal Freedom ±45° 
m^ Brighter than +3.5, Real Stars 


Initial Att. 
Covariance 
T1/T2/T3 
(arc-sec) 

Initial Gyro Bias 
Covariance 
w 

(deg/hr) 

Sensor Noise 
SN 

(arc-sec) 

Gyro Noise 
SV 

Gyro Noise 
SU 

10/10/10 

.05 

3 

.115 

3.33x10“^ 





(2 sic in 

(.002 deg/hr/hr) 





5 minutes) 


lO/lOy 

'10 

.05 

3 

.05 

3.2x10'^ 



.05 

3 

.05 

1.6x10“^ 



.05 

3 ' 

.015 

3.2x10"^ 



.05 

3 

.015 

1.6x10“^ 

360/3C 

50/360 

.05 

3 

.05 

1.6x10"^ 

10/10/10 

0.1 

3 

.05 

1.6cl0"^ 

360/360/360 

0.1 

3 

.05 

1.6x10"^ 


T1,T2^T3 = 10.0 10.0 iP,0 

M = 5,00E-02 

SN,SV,SU = 3.00E+00 1.15E-01 3.33E-05 

H = 19300.0 


SI S2 S3 PS RS Tt T2 T3 PT «T 


0.0 

394.3 

798.6 

0. 0 
22.2 
11, 0 

0. 0 
22,2 
26. 7 

2 

3 

1434.3 

8.9 

22. 1 

1 

1478.8 

3.1 

22.9 

1 

1479.8 

2.2 

2. 0 




5849.0 ro 
6787.7 


9.6 

10.1 

10.1 



.6 

.3 
.0 

10.1 

11.2 

10.6 


3 

0 

1 


5.6 

5.0 

20.4 
12.9 

13.5 
13.2 


4.4 
3.2 

12.6 

9.5 

12.1 

3.9 




















30101.5 

2.7 

2.5 

6.8 

3,7 

7.8 

2.0 

1.4 

2.4 

2.5 

3.5 


31081,5 

4 . 8 

h * 5 

5.1 

6.6 

8 ♦ 3 

2.6 

2.9 

3.9 

3.8 

5.5 


31092.3 

2.6 

2. 9 

4.0 

3.9 

5.6 

2.0 
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